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Abstract. Niche breadth and niche overlap measurements have been widely used together to assess resource dynamics in 
biological communities. However, the estimations of these measurements are susceptible to bias due to sample size. Al- 
ternatives such as standardization of niche breadth have been implemented to try and mitigate this bias, but standardiza- 
tion does not solve the underlying problem. In anuran trophic ecology, sampling is usually incomplete, which constrains 
comparisons. A solution to this problem could lie in implementing an analysis using H1Lv’s series (also called HILL num- 
bers) for comparisons under the same sample coverage; however, this approach has not yet been tested in anuran trophic 
ecology research. The objectives of this work were to evaluate the bias resulting from sample size for the most frequently 
used measurements of niche breadth and niche overlap in anuran trophic ecology vis-a-vis overlap measurements derived 
from HI11’s series and to provide a least biased protocol for anuran trophic ecology evaluations. We built data matrices 
with similar features to those expected for anuran assemblages and quantified the bias of each measurement for different 
sample sizes. We found that HILL’s series measurements were less biased and more informative than traditional ones. We 
provide an analysis guideline based on HiLv’s series that facilitates direct comparisons between predator species regarding 
their consumed and shared prey communities. 
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Introduction 


Niche breadth and niche overlap have been used as com- 
plementary measurements in descriptions of resource use 
in biological communities (KREBS 1989, SOLE & RODDER 
2010). However, both of these have limitations, including 
the scale at which resources are measured, the spatial dis- 
tribution of samples, the way that resources are categorized 
(COLWELL & FUTUYMA 1971, DE CACERES et al. 2011), and 
the sample size effect (RICKLEFS & LAU 1980, SMITH & ZA- 
RET 1982, Kovács & TÖRÖK 1997). The latter is especially 
important because it may change the biological interpreta- 
tion of the niche breadth and overlap measurements (e.g., 
underestimated niche breadth or overlaps in species as- 
semblages due to low sample sizes) (KREBS 1989, KovAcs 
& TOROK 1997). Although the bias resulting from sample 


size niche breadth measurements was not yet evaluated at 
the time, KovAcs & TOROK (1997) suggested that for niche 
breadth comparisons, regardless of the measurement used, 
a minimum number of samples for each species should 
be set for their biological interpretation to be meaningful. 
Based on simulations, SMITH & ZARET (1982) then evaluat- 
ed the sample size effects on seven niche overlap measure- 
ments, for which the MorisiTa and RENKONEN measure- 
ments were the least and the most biased, respectively, in 
the niche overlap estimation at minimum sample size (~40 
samples). These authors attributed this bias to the weight 
lent by each measurement to the resource categories with- 
in the sample, and to the evenness of these categories in 
the sample. It is important to note that SMITH & ZARET 
(1982) excluded some of the overlap measurements most 
frequently used in ecology research, such as PIANKA’s sym- 
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metric index, which is the most commonly used measure- 
ment in trophic ecology research of amphibian and reptile 
communities, and for which the estimated bias due to vari- 
ations in sample size has not been evaluated (TOFT 1980, 
KREBS 1989, SOLE & RODDER 2010). 

Trophic interactions in anuran communities are useful 
ecological models for assessing bias in niche breadth and 
overlap measurements, because compared to other verte- 
brate groups, anuran diets generally show (i) a compara- 
tively low number of prey animalsconsumed per individu- 
al (X ~9.85); (ii) a relatively low number of prey categories 
per species (X ~9.75); and (iii) that the minimum number 
of samples needed for a complete sampling is around 80- 
115. These sample sizes are often not achieved in short-term 
samplings because most species’ activity periods are restrict- 
ed to specific seasons and are closely related to climatic vari- 
ations (TOFT 1980, ZUG 1993, KOVÁCS & TÖRÖK 1997, PARM- 
ELEE 1999, MANEYRO et al. 2004, VITT & CALDWELL 2014). 

Niche breadth measurements are conceptual applica- 
tions of alpha diversity measurements (KREBS 1989). The 
difference between alpha diversity and niche breadth meas- 
urements is that the latter are not used to quantify spe- 
cies diversity (as are alpha diversity measurements), but 
the richness of resource types used by a taxon (COLWELL 
& FUTUYMA 1971, KREBS, 1989). On the other hand, niche 
overlap measurements are a beta diversity measurement 
application (KREBS, 1989). The difference is that overlap 
measurements are not used to quantify distances in species 
composition and abundance between two or more commu- 
nities, but quantify the proportion of resources shared by 
two or more taxa (COLWELL & FUTUYMA 1971, KREBS 1989). 

The most frequently used niche breadth measurements 
for anuran diet analyses are the SHANNON and LEVINS 
measurements (LEVINS, 1968, COLWELL & FUTUYMA 1971, 
SOLE & RODDER 2010). LEVINS’ measurement is equal to the 
inverse of Srmpson’s index, which, like SHANNON’s meas- 
urement, is an index used to evaluate community struc- 
ture. However, applying these indices (based on relative 
abundances) in trophic ecology research is fraught with 
the same problem as assessing the diversity of a biological 
community: they lack biologically interpretable units. For 
example, SHANNON index units depend on the logarithm 
used. Thus, the magnitude of the SHANNON index does not 
explicitly reveal how the community is structured accord- 
ing to the distribution of relative abundances across species 
or resource items (e.g., prey) (see MACARTHUR 1965). The 
absence of biological units (e.g., species or operational tax- 
onomic units, OTUs) also limits unbiased comparisons be- 
tween units of analysis (see below). HILL (1973) addressed 
the cautioning remarks by MACARTHUR (1965) and pro- 
posed that the diversity or structure of a community (i.e., 
how relative abundance is distributed among species) be 
assessed from a mathematical series that has since become 
known as Hi1Lt’s numbers or the H11r’s series. The HIL1’s 
numbers formula is presented in detail in the Materials 
and methods section of this paper, but it may be helpful to 
highlight that: i) H1Lu’s numbers allow to estimate commu- 
nity diversity in biologically interpretable units (i.e., the ef- 
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fective number of species); ii) HILL’s main numbers or di- 
versity orders are N = o, which estimates species richness, 
N = 1 is the effective number of equally common species, 
equivalent to the exponential of SHANNON’s index, and N = 
2 estimates the effective number of abundant or dominant 
species, equivalent to the reciprocal of Srmpson’s index; 
and iii) that they offer a unified mathematical framework 
for the comparison of diversity among communities (HILL 
1973, JOST 2006). 

The most used overlap measurement in anuran trophic 
ecology research is PIANKA’s (1973) symmetric measure- 
ment, which is mathematically similar to MacARTHUR’s 
and LEvINs’s overlap measurement, except that the latter 
is asymmetrical (PIANKA 1974, KREBS 1989, SOLE & RÖD- 
DER 2010). Other measurements less widely employed for 
trophic niche analysis are RENKONEN’s percentage overlap 
or Schoener’s measurement (KREBS 1989), and MORISI- 
TA’s similarity measurement, which has been suggested to 
be the least biased by sample size (SMITH & ZARET 1982). 
Additionally, since the niche breadth measurements most 
frequently used in anuran dietary analysis are derivatives 
of the orders 1 and 2 of HILr’s series, the species turnover 
measurements for the diversity orders 1 and 2 (Jost 2007) 
could be used. 

A recurrent aim in trophic ecology studies of anuran 
communities is the comparison between species (e.g., 
TOFT 1980, PARMELEE 1999, SOLE & RODDER 2010). How- 
ever, niche breadth measurements cannot be compared di- 
rectly, because their value ranges depend on the number 
of each species’ resources (KREBS 1989). To address this 
problem of comparability, many authors have proposed the 
standardization of niche breadth measurements, adjusting 
their value range between o and 1 (KREBS 1989, MANEY- 
RO et al. 2004, LIMA, RODDER & SOLE 2010, SOLE & RÖD- 
DER 2010, FORTI et al. 2011, MARf{N-MARTINEZ et al. 2019). 
However, since the most commonly used measurements of 
niche breadth (i.e., SHANNON and Levins) can be obtained 
from HILv’s series, comparisons can (and should) be made 
considering the concept of sample coverage (CHAO & JOST 
2012). 

Sample coverage can be understood as “the actual rela- 
tive abundance of the community represented by the pool 
of species in a sample” and can be used as a criterion for the 
comparison of samples at the same level of completeness 
through an interpolation/extrapolation protocol (CHAO 
& JOST 2012, CHAO et al. 2014). Furthermore, when com- 
paring measurements at the same level of sample cover- 
age, the sample size bias is minimized by default (CHAo & 
Jost 2012, ENGEL et al. 2021), allowing direct comparison 
of niche breadth measurements between species independ- 
ent of species abundance. However, this concept has not 
yet been used in anuran trophic ecology studies (or in any 
other biological group, to our knowledge). The sample cov- 
erage approach should be applied to comparisons of beta 
diversity, too (overlap in the case of anuran diets). Howev- 
er, how exactly this concept should be applied in this area 
is not fully resolved, and its application is still under devel- 
opment (ENGEL et al. 2021). 
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Our general objective was to evaluate the bias due to 
changes in sample size in the measurements most frequent- 
ly used in dietary assessments of anurans and compare it 
with values obtained from the HIL1’s series framework. 
Considering that sample size effects have not yet been eval- 
uated for niche breadth measurements but been limited to 
a few niche overlap measurements, and given that achiev- 
ing a minimally informative sample number for the study 
of anuran trophic ecology is a central challenge, we ad- 
dress the following questions: (i) how does the ecological 
interpretation differ in the analysis of prey communities of 
frog species using traditional niche breadth measurements 
from the HiL1’s series approach’; (ii) how do the ecologi- 
cal interpretations differ in the analysis of the trophic niche 
overlap of a frog assemblage using a traditional overlap and 
HILL’s series approach?; (iii) how does the bias in the esti- 
mation of trophic niche breadth and overlap for simulated 
anuran species change with variations in sample size?; and 
(iv) which are the least biased breadth and overlap meas- 
urements for dietary descriptions of anurans? Therefore, 
we here provide a suitable dietary analysis protocol for an- 
urans considering their trophic characteristics and sensi- 
tivity to sample size. 


Materials and methods 
Comparison of approaches 


To compare traditional and HILL’s series measurements 
to evaluate niche breadth and overlap, we used diet data 
obtained for an assemblage of four frog species inhabiting 
avocado orchards in central-western Mexico (Eleuthero- 
dactylus angustidigitorum, E. nitidus, Dryophytes arenicolor 
and Lithobates forreri). Details of these diets can be found 
in the Table S1. We considered as traditional niche breadth 
and niche overlap measurements those compiled by SoLÉ 
& RODDER (2010). For niche breadth we used two measure- 
ments, the standardized SHANNON’s measurement, which 
is mathematically equal to PrELou’s measurement of even- 
ness: 


[= —Xpilog pi (1) 
logn 


where 'J": standardized SHANNON’s measurement; p: 
proportional abundance of each prey category; and n: 
number of prey categories; and the standardized Levins’s 
niche breadth measurement, which is equivalent to the 
standardized Simpson’s index: 


p- a" (2) 


n-1 


where B He standardized LEVINS’ measurement, p: rela- 
tive abundance of each prey or category; and n: number of 
prey categories. 


As traditional niche overlap measurement, we used P1- 
ANKa’s symmetric index which can be calculated as: 


Oj L È PijPik 
EP pA? pi 6) 
where O, is the niche overlap and p, and p, represent 


the i™ resource proportion used by 'j' and 'k' species. 

Niche breadth measurements were compared directly. 
For their calculation we used the R environment package 
'abdiv' (BITTINGER 2020). The PIANKA’s symmetric index 
was calculated using the R environment package 'spaa' 
(ZHANG 2016). 

For description of prey communities based on HILL'S 
series, we applied the Jost (2006) mathematical notation: 


D= (Chip = (Mw 


where p,is the i prey proportional abundance, S is the 
prey richness, and q is the diversity order (for orders differ- 
ent from 1). When q = 1, diversity should be calculated as: 


1D = exp(D}_, pj log p;) (5) 


and it is equivalent to the exponential SHANNON’s meas- 
urement. 

To compare prey diversity among anuran species, we 
used o, 1 and 2 diversity orders (i.e., °D, ‘D and °D), and the 
comparison was carried out with equal sample coverage: 


¢, = 1-2[-S | (6) 
n L(n-1)f,+2f 

where Cn is the sample coverage, f, is “singletons” (prey 
with abundance equal to 1), f, is “doubletons” (prey with 
abundance equal to 2), and n is the total number of prey- 
registered individuals. Sample coverage allows the estima- 
tion of the proportion of the community represented by 
the prey in the sample; using the same sample coverage 
level allows a direct comparison of diversity order values 
(CHAO & JosT 2012). To evaluate differences in the estima- 
tions, we used comparisons on the 95% confidence inter- 
val overlaps (CUMMING, FIDLER & VAUX 2007), estimated 
via bootstrap (10,000 repetitions) (LEGENDRE & LEGENDRE 
2012). The 4D + 95% CI per frog species was estimated us- 
ing the NEXT" package of R (Hsu et al. 2016). 

For overlap measurements, we used 1 - 49D, in orders o, 
1and 2 (i.e., “Dp 'D, and Dy modified from Jost (2007): 


1/Q- vO 
Dg= 1—[(%/ Ma) |= Ba? Gy 
where A is the order q alpha diversity measurement 
(local diversity, in this case single-species prey diversity), 
IÀ is the order q gamma diversity measurement (landscape 
diversity, in this case all-species prey diversity), and q is 
the diversity order (for orders different from 1). When q = 
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1, beta diversity should be calculated as the q limit to 1 (see 
Jost 2007 for a wider explanation). For the 9D, calculation, 
we used the ‘vegetarian’ package of the R environment 
(CHARNEY & RECORD 2012). 

To evaluate differences in estimations of overlaps be- 
tween species pairs, we conducted visual comparisons 
based on 95% confidence intervals (CUMMING, FIDLER & 
VAUX 2007), estimated via bootstrap (10,000 repetitions) 
(LEGENDRE & LEGENDRE 2012). 


Selecting niche breadth and overlap measurements 


The niche breadth measurements evaluated in this work 
were both of the orders 'D (equivalent to SHANNON’s ex- 
ponential measurement) and °D (equivalent to LEvINs’ 
niche breadth measurement and the inverse of Simpson’s 
index) proposed by Hitt (1973). For measuring overlaps, 
we used PIANKA’s symmetric measurement, overlap per- 
centage (RENKONEN’s measurement), MorisitT«’s (similar- 
ity) overlap measurement (KREBS 1989), and 1 - turnover 
measurements for diversity orders 1 and 2 of the HILL’s se- 
ries framework (Jost 2007). For all measurements, species 
x prey category sum vectors were considered, which were 
constructed from anuran individuals x prey category input 
matrices. 


Input matrices 


For niche breadth, we considered hypothetical species with 
two known trophic niche breadth values for each measure- 
ment, an intermediate ('D: 6.42; °D: 5.04) and a maximum 
value, equal to prey item richness ('D: 10; *D: 10). We con- 
sidered species with 10 prey categories and a mean prey 
number of 10 individuals. For this purpose, we used the 
'rpois' function of the R environment (R Core Team 2021) 
with lambda = 10 and n = 10 to generate a vector of 10 ran- 
dom numbers (prey categories). From the generated vector 
we constructed abundance matrices of 200 anuran individ- 
uals (maximum sample size). 

For overlap measurements, two hypothetical species 
with 15 prey categories in total (10 for each species) were 
considered in two scenarios: (i) five shared prey catego- 
ries (overlap ~50%); (ii) all prey categories shared (overlap 
~100%). We considered a mean prey number of 10 indi- 
viduals and used the 'rpois' function (R Core Team 2021) 
to generate two vectors with a 15 random numbers vector 
(prey categories) for each frog species. From the generated 
vectors we constructed two frog species abundance matri- 
ces of 200 individuals each (maximum sample size per spe- 
cies considered here). 


Simulation procedure 


The following procedure was used both for the amplitude 
and niche overlap simulation protocols. Sum vectors were 
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generated by randomly choosing anuran individuals from 
the input matrices, with sample sizes of 5, 10, 20, 40, 80, 160 
and 200 individuals. We calculated 5,000 times the meas- 
urement for each sample size using random sum vectors. 
Then, we estimated: (i) the mean of the estimation for each 
sample size; (ii) the standard deviation for each estima- 
tion for each sample size; (iii) the 95% confidence intervals 
for each sample size; and (iv) the bias of the estimation for 
each sample size, calculated as: 


observed value (8) 


bias = 1 - 
expected value 


where the ‘observed value’ was the value of the mean for 
each sample size, and the ‘expected value’ was the known 
value of the sum vector initially generated. This bias calcu- 
lation takes values between o and 1 and can be transformed 
to a percentage. All analyses and the functions used for our 
simulations were built in the R environment version 4.1.2 
(R Core Team 2021), using the ‘vegetarian’ (CHARNEY & 
RECORD 2012) and 'spaa’ (ZHANG 2016) packages. 


Results 


The results when comparing the trophic niche breadth of 
the frog assemblage using the traditional and the H11t’s 
series approaches were different (Fig. 1). The traditional 
approach detected no differences between species’ niche 
breadths (Figs 1A-B). The Hixt’s series-derived approach 
detected differences between species’ trophic niches at °D, 
'D and ’D (Figs 1C-E). Regarding trophic niche overlap 
measurements, there were no significant differences be- 
tween species pairs using the PIANKA’s symmetric index 
(Fig. 2A). Likewise, there was no significant difference be- 
tween species pairs using `D, and *D,, while °D, detected 
statistically significant differences between Dryophytes 
arenicolor and Lithobates forreri and the other species pairs 
(Figs 2B-D). 

All niche breadth and overlap measurements evaluated 
in this work exhibited biases with respect to sample size 
(Fig 3. 3, 4, Table S2). As for niche breadth, with the 'D (i.e. 
SHANNON’s exponential measurement), the maximum bias 
for the expected value of 6.42 was 80.70% (Fig. 5A) and 
for the expected value of 10 the maximum bias was 77.07% 
(Fig. 5B). With °D (Levins niche breadth measurement), 
for the expected value of 5.04 the bias was 82.53% (Fig. 5C), 
and for the expected value of 10 the maximum bias was 
78.52% (Fig. 5D). 

Likewise, the sample size at which the estimates were 
stabilized (i.e., where the bias of the expected value is less 
than 10%) was 80 samples for expected value 10 and 160 
samples for expected value 6.41 of the 'D measurement, 
and 160 samples for both expected values 5.04 and 10 of 
the 7D measurement ( 5). 

For the niche overlap measurements, the maximum bias 
occurred at low sample size values (5-10 samples, Figs 4, 
6), and ranged between 58 and 82%. In decreasing bias or- 
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der, the less biased measurements were: °D, (58.57% for 
0.5 expected overlap, and 60.55% for 1 expected overlap; 
Figs 6I-J), 'D, (64.96% for 1 expected overlap, and 65.82% 
for 0.5 expected overlap; Figs 6G-H), PIANKA’s symmetric 
measurement (71.97% for 1 expected overlap, and 76.21% 
for 0.5 expected overlap; Figs 6E-F), MortsiTa’s measure- 
ment (72.94% for 0.5 expected overlap, and 82.27% for 1 ex- 
pected overlap; Figs 6C-D) and RENKONEN’s overlap per- 
centage (76.47% for 0.5 expected overlap, and 81.25% for 1 
expected overlap; Figs 6A-B). 

The overlap measurements that stabilized faster (first 
reaching bias values of less than 10%) with increasing sam- 
ple size were: 'D, (40 samples for 1 expected overlap, and 80 


samples for 0.5 expected overlap), 7D, and PIANKa’s sym- 
metric measurement (80 samples, for both 0.5 and 1 ex- 
pected overlap), MorisiTa’s measurement (80 samples for 
1 expected overlap, and 160 samples for 0.5 expected over- 
lap) and RENKONEN’s overlap percentage (160 samples for 
both 0.5 and 1 expected overlap; Fig. 6). 


Discussion 
Hirr’ series for trophic ecology analysis of anuran assem- 


blages allow a more comprehensive and realistic biologi- 
cal view of the observed patterns than do those obtained 
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Figure 1. Evaluation of the prey community consumed by an anuran assemblage in avocado orchards in central-western México 
(mean + 95 % CI). (A) SHANNON’s standardized measurement of niche breadth; (B) Levins’ standardized measurement of niche 
breadth; (C) °D of Hixt’s series, equivalent to prey richness; (D) 'D from HILL’s series, equivalent to the effective number of equally 
common prey items; (E) 7D of H1L1’s series, equivalent to the effective number of dominant prey. EA: Eleutherodactylus angustidigito- 
rum, EN: Eleutherodactylus nitidus, DA: Dryophytes arenicolor, LF: Lithobates forreri. 
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Figure 2. Overlap evaluations in a frog assemblage in avocado orchards in central-western México through distinct methods (mean + 
95 % CI). (A) PIANKa’s symmetric measurement; (B) "D,, equivalent to overlap by prey richness; (C) 'D,, equivalent to overlap by 
equally common prey; (D) *D,, equivalent to overlap by dominant prey. EA: Eleutherodactylus angustidigitorum, EN: Eleutherodactylus 
nitidus, DA: Dryophytes arenicolor, LF: Lithobates forreri. 
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Figure 3. Variations in the trophic niche breadth estimation with respect to sample size. Left: Expected values 6.41 above and 5.06 
below. Right: Expected value for both cases 10. (A-B) 'D (i.e., SHANNON’s exponential); (C-D) °D (i.e., Levins’ niche breadth meas- 
urement, cf. SImMpson’s inverse index). Grey line: Expected value for each case. 
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through the use of traditional measurements of niche 
breadth and overlap. Our results suggest that all niche 
breadth and niche overlap measurements evaluated were 
biased at small sample sizes (5-10 samples) and that bias 
decreased when sample size increased. For niche breadth 
measurements, the percentage of bias difference was less 
than 3%. For overlap, the bias percentage differed between 
measurements, with 'D, and °D, being the least biased ones 
and their difference to the next least biased measurement 
(PIANKA’s symmetric measurement) was approximately 
7%, and about 10% to the most biased overlap measured 
(RENKONEN measurement). Regarding estimation stabili- 
zation (a sample size where the bias is less than 10%), the 
earliest niche breadth measurement to stabilize was 'D, 
while for the overlap one it was 'D,. 


A B 
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Figure 4. Mean overlap estimations across different sample sizes 
for different measurements. Left: Estimate when half of the com- 
position is shared (overlap ~ 50%). Right: Estimation when all 
of the composition is shared (overlap~100%). (A-B) Percentage 
overlap (RENKONEN, 1938); (C-D) Morisita (1959); (E-F) P1- 
ANKA’s symmetrical measurement (PIANKA, 1973); (G-H) 'D, 
(Jost 2007); (I-J) *D, (Jost 2007). Grey line: Expected value for 
each case. 


The interpretations of trophic analyses of anurans from 
the traditional approach were different from those using 
the HiL1’s series approach. According to the observations 
with the traditional approach, there were no significant 
differences between the trophic niche amplitudes in any 
of the species in the anuran community, while in the ap- 
proach using the HILv’s series and the comparison based 
on sampling coverage, differences in the prey communi- 
ty consumed by each species were observed in three ways: 
(i) in the prey richness consumed by the species of the as- 
semblage, forming three different groups; (ii) three differ- 
entiated groups of species were formed according to their 
equally common prey; and (iii) two differentiated groups 
of species were formed according to the effective number 
of their dominant prey. 

Kress (1989) and SOLE & RODDER (2010) suggested em- 
ploying a standardized trophic niche breadth (mathemati- 
cally adapted so that their value domain stays between o 
and 1) to compare species. However, standardization limits 
the actual variation observed and may cause misinterpre- 
tation when comparing the niche breadth between two or 
more species. Trophic niche breadth measurements are al- 
pha diversity measurements (i.e., orders 1 and 2 of HILL'S 
series), so that CHAo et al’s (2014) suggestion for alpha di- 
versity comparisons is applicable to them: the communi- 
ties (species in this case) should be compared at the same 
percentage of sample coverage to avoid biological misin- 
terpretations due to sample size effects. The use of sample 
coverage in anuran diet analysis permits the comparison of 
anuran diet diversity (i.e., niche breadth) even if their sam- 
ple sizes are different, if they are at the same level of sample 
coverage (KovACs & TÖRÖK 1997, CHAO et al. 2014). 

The bias between overlap measurements was variable. 
The 'D, and °D, estimators were less biased and required 
fewer samples to reach less than 10% bias. On the other 
hand, the MorisiTa’s measurement and the RENKONEN’S 
overlap percentage were the most biased estimators with 
respect to sample size. These results contrast with SMITH 
& ZARET (1982), who found that the MorisitTa’s meas- 
urement was the less biased measurement in the estima- 
tion of overlap with respect to sample size. This may be 
because SMITH & ZARET (1982) used a minimum sample 
size of 50 samples, while our study used a minimum sam- 
ple size of 5 samples, and we found that MorisiTa’s meas- 
urement bias decays by between 40-60% when sample size 
increases from 5 to 50 samples. SMITH & ZARET (1982) and 
KREBS (1989) suggested that the most biased overlap meas- 
urement with respect to sample size is RENKONEN’s over- 
lap percentage. This result was consistent with our findings 
in that RENKONEN’s overlap percentages had the highest- 
percentage biases of all overlap measurements with respect 
to sample size and was the measurement that required the 
largest number of samples to stabilize (i.e., to reach a bias 
value of less than 10%). 

PIANKA’S (1974) statement that the use of distinct over- 
lap measurements is “somewhat arbitrary since similar 
qualitative results are obtained with a wide variety of indi- 
ces” can be considered correct, because all measurements 
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underestimate overlaps at small sample sizes (5-10 sam- 
ples), but this underestimation (bias) was not the same 
for all overlap measurements. Our data suggest that both 
Pranxka’s and Morisita’s overlap measurements should 
not be used for sample sizes smaller than 80 individuals 
and RENKONEN’s overlap measurement should not be used 
with sample sizes smaller than 160 individuals. For sam- 
ple sizes smaller than 80, the most recommendable meas- 
urements (due to their lower biases) are 'D, and D, (see 
above). 

We found differences in the ecological interpretations 
derived from traditional analysis approaches to the de- 
scription of prey communities in anuran assemblages and 
in their overlap analyses. In general, the HIL1’s series ap- 
proach facilitated visualization and comparison in a man- 
ner that traditional measurements did not (i.e., richness, 
equally common prey and dominant prey), so that inter- 
pretation and the description of anuran trophic interac- 
tions become more thorough. Additionally, due to the im- 
plementation of the comparison by sampling coverage, a 
sample-size bias in comparisons is avoided (CHAO & JOST 
2012), allowing prey communities to be compared directly 
and without the need to use standardizations that distort 
observed patterns. 

The 'D and °D measurement biases were similar. How- 
ever, since these are equivalent to orders 1 and 2 of the 
Hitv’s series, their use should not be mutually exclusive, 
but rather complementary, because each one reflects a dif- 
ferent portion of diversity (trophic niche) (HILL 1973, Jost 


A 
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2006). Regarding the interpretation of each measurement 
(i.e., order of diversity) in an anuran dietary context, ‘D 
would reflect the equally abundant dietary components, 
without bias for rare or abundant prey, and ?D reflects the 
most abundant (i.e., dominant) dietary components in the 
prey community that are consumed by the anuran species 
(JOST 2018, CULTID-MEDINA & ESCOBAR 2019). 

Due to sample size bias results found in both niche 
breadth and niche overlap measurements, and because 
HILL’s series approach provides a more adequate trophic 
niche quantification and biological interpretation, we pro- 
pose the following dietary analysis guidelines for anurans: 

1. Prey diversity consumed by anuran species (trophic 
niche) should be described using the diversity measure- 
ments °D: Consumed prey richness, ‘D: Number of equally 
common prey (equivalent to SHANNON exponential meas- 
urement), and *D: Effective number of dominant prey 
items (equivalent to LEvins niche breadth measurement). 

2. Comparison of any of the diversity orders (i.e., °D, ‘D 
and °D) should be carried out when species are at the same 
level of sample coverage, or by interpolating or extrapo- 
lating their values to ensure that the comparison is at the 
same level (see CHao et al. 2014 for an explanation of the 
interpolation/extrapolation protocol). 

3. Anuran dietary similarity (overlap) should be com- 
pared using as measurements the inverse Jost (2007) turn- 
over (i.e., 1-2D,) of orders o, 1 and 2 as follows: °D; overlap 
by prey richness, 'Dy overlap by equally common prey, and 
*D,;: overlap by dominant prey. 


Sample size 


Figure 5. Bias in the estimation of trophic niche breadth as a result of sample size. Bias is estimated as abs (1 - estimated value/expected 
value). (A-B) 1D (i.e, SHANNON’s exponential); (C-D) °D (i.e., Levins’ trophic niche breadth measurement, i.e,. inverse SIMPSON 
index). Grey line: Expected value for each case, dotted line: 10% bias. 
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4. Although overlap measurements in anuran diets 
should be compared when the sampling effort (sample 
coverage) is the same for all species (as beta diversity com- 
parisons), the manner in which sample coverage should 
be applied to beta comparisons is not fully resolved (EN- 
GEL et al. 2021), and its resolution is beyond the scope of 
this paper. 

5. Niche breadth and overlap measurements should be 
accompanied by complementary and more descriptive 
ones such as the index of relative importance (IRI) (PINKAS 
et al. 1971; HART et al. 2002), to arrive at more accurate 
and complete interpretations of the trophic dynamics of 
anuran communities than those obtained by measuring 
trophic niche breadth and overlap only. 


150 200 0 50 100 


150 200 


100 150 200 0 50 100 150 200 


Sample size 


Figure 6. Bias in overlap estimation for different measure- 
ments at different sample sizes. The bias is estimated as abs 
(1 - estimated value/expected value). Left: Estimation when 
half of the composition is shared (overlap ~50%), right: Esti- 
mation when all of the composition is shared (overlap ~100%). 
(A-B) Percentage overlap (RENKONEN 1938); (C-D) MORISITA 
(1959); (E-F) P1aAnKa’s symmetrical measurement (PIANKA 
1973); (G-H) "Ds (Jost 2007); (I-J) D; (JosT 2007). Dotted 
line: 10% bias. 
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